Pattern inspection apparatus, image alignment method, displacement amount estimation method, and computer-readable recording medium with program recorded thereon

ABSTRACT

A pattern inspection apparatus includes a first unit configured to acquire an optical image of pattern, a second unit configured to generate a reference image to be compared, a third unit configured to calculate elements of a normal matrix for a least-squares method for calculating a displacement amount displaced from a preliminary alignment position, a forth unit configured to estimate a type of the reference image pattern, by using some of the elements of the normal matrix, a fifth unit configured to calculate the displacement amount based on the least-squares method, by using a normal matrix obtained by deleting predetermined elements depending upon the type of the pattern, a sixth unit configured to correct an alignment position between the optical image and the reference image to a position displaced by the displacement amount, and a seventh unit configured to compare the optical image and the reference image.

CROSS-REFERENCE TO RELATED APPLICATION

This application is based upon and claims the benefit of priority fromthe prior Japanese Patent Application No. 2006-217672 filed on Aug. 10,2006 in Japan, the entire contents of which are incorporated herein byreference.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to pattern inspection apparatus, an imagealignment method, a displacement estimation method or a computer programto apply the method. For example, it is related with a technique forpattern inspection which inspects a pattern defect of a template used inmanufacturing semiconductor devices, and an apparatus which inspects adefect of a micro-pattern of a photomask, a wafer, or a liquid crystalsubstrate used in manufacturing a semiconductor device or a liquidcrystal display (LCD).

2. Description of Related Art

In recent years, with increasing integration of large volume oflarge-scale integrated (LSI) circuits, circuit line-widths required forsemiconductor devices are becoming narrower and narrower. Thesesemiconductor devices are manufactured by transferring a photo-maskpattern onto a wafer by means of a reduced-magnification-projectionexposure apparatus (a stepper) while using a template pattern (called amask or a reticle, and will be called a mask hereinafter) on whichcircuit patterns are written. Therefore, a pattern writing apparatuswhich can write micro-circuits onto a mask is deployed in writing a maskfor transferring micro-circuit patterns onto a wafer. A pattern circuitmay be directly written onto a wafer by using the pattern direct-writingapparatus. In addition to a writing apparatus using electron beams, alaser beam pattern writing apparatus which uses laser beams to write apattern is also developed.

Improvement in yield is crucial in manufacturing an LSI which requires alot of cost. However, as one-gigabit DRAM (Random Access Memory), theprecision of a pattern has been changing from sub-microns to nanometers.One of major factors which decrease the yield is pattern defects of amask pattern used in exposing and transferring a micro-pattern onto asemiconductor wafer by a photolithography technique. In recent years,with miniaturization of an LSI pattern written on a semiconductor wafer,dimensions which have to be detected as a pattern defect are becomingextremely small. Therefore, a pattern inspection apparatus whichinspects defects of a transfer mask used in manufacturing an LSI needsto be highly precise.

On the other hand, with development of multimedia, the size of a liquidcrystal substrate of an LCD (Liquid Crystal Display) is becoming large:500 mm×600 mm or more, and miniaturization of a pattern of a thin filmtransistor (TFT) or the like formed on a liquid crystal substrate isadvancing. Therefore, it is required that a considerably small patterndefect should be inspected in a large area. For this reason, developmentof pattern inspection apparatus which efficiently inspects a defect of apattern of a large-area LCD and a photomask used in manufacturing thelarge-area LCD in a short time is urgently required.

As to a conventional pattern inspection apparatus, it is well-known thatan inspection is performed by comparing an optical image captured byphotographing a pattern written on pattern, such as a lithography mask,at a predetermined magnification by using a magnifying optical systemwith design data or an optical image captured by photographing the samepattern on the target workpiece (see, e.g., Published UnexaminedJapanese Patent Application No. 08-76359 (JP-A-08-76359)).

For example, the following is known as pattern inspection methods:die-to-die inspection which compares optical image data obtained bycapturing the same patterns at different positions on the same mask, anddie to database inspection which inputs drawing data (design patterndata) obtained by converting CAD data into appropriate format to beinputted by a drawing apparatus when drawing a pattern on a mask, intoan inspection apparatus, generates design image data (reference image)based on the inputted drawing data, and compares the generated designimage data with an optical image serving as measurement data obtained bycapturing an image of the pattern. In these inspecting methods of theinspection apparatus, pattern is placed on a stage to be scanned by aflux of light when the stage moves to perform inspection. The targetworkpiece is irradiated with flux of light from a light source and anillumination optical system. Light transmitted through the targetworkpiece or reflected by the target workpiece is focused onto a sensorthrough an optical system. The image captured by the sensor istransmitted to a comparing circuit as measurement data. In the comparingcircuit, after alignment of the images, the measurement data is comparedwith reference data based on appropriate algorithm. When the measurementdata is different from the reference data, it is estimated there to be apattern defect.

Herein, the reference image and the optical image are compared in eacharea of a predetermined size. Highly precise alignment between thereference image and the optical image is required for performing thiscomparison. A technique for calculating displacement amount or“deviation” between a reference image and an optical image by use of aleast-squares method is disclosed in a reference (for example, refer toJP-A-11-153550). Further, an interpolation method for interpolatingimage data to be obtained by use of neighboring 4-point or 16-pointimage data is described in a reference (for example, refer to ImageAnalysis Handbook, pp. 442 to 443, University of Tokyo Press, firstedition issued on Jan. 17, 1991).

With a miniaturization of a pattern, there is a demand for a furtherprecision of alignment required for detecting micro-defects. The pointherein is to correct only systematic error factors such as a stageplacement error, a speed error or a magnification error, but not tocorrect inconsistent portions that occur locally and randomly such asdefects, if possible.

As described above, highly precise alignment between a reference imageand an optical image is required for performing comparison. However,with a miniaturization of a pattern, it has become difficult to detectrelative displacement between the reference image and the optical imagein high precision.

BRIEF SUMMARY OF THE INVENTION

It is an object of the present invention to provide a method and anapparatus which perform highly precise positional alignment between areference image and an optical image.

In accordance with one aspect of the present invention, a patterninspection apparatus includes an optical image acquiring unit configuredto acquire an optical image of pattern to be inspected on which apattern is formed, a reference image generating unit configured togenerate a reference image to be compared with the optical image, anormalal matrix element calculating unit configured to calculate aplurality of elements of a normalal matrix for a least-squares methodfor calculating a displacement amount displaced from a preliminaryalignment position between the optical image and the reference image, apattern type estimating unit configured to estimate a type of a patternindicated by the reference image, by using some of the plurality ofelements of the normalal matrix, a displacement calculating unitconfigured to calculate the displacement amount displaced from thepreliminary alignment position based on the least-squares method, byusing a normalal matrix obtained by deleting predetermined elements fromthe plurality of elements of the normalal matrix depending upon the typeof the pattern which has been estimated, a position correcting unitconfigured to correct an alignment position between the optical imageand the reference image to a position displaced from the preliminaryalignment position by the displacement amount, and a comparing unitconfigured to compare the optical image and the reference image whosealignment position has been corrected.

In accordance with another aspect of the present invention, a patterninspection apparatus includes an optical image acquiring unit configuredto acquire an optical image of pattern to be inspected on which apattern is formed, a reference image generating unit configured togenerate a reference image to be compared with the optical image, afirst SSD (Sum of Squared Difference) calculating unit configured tocalculate a first displacement amount from a preliminary alignmentposition between the optical image and the reference image to a firstposition where an SSD between a pixel value of the optical image and apixel value of the reference image is minimized, a normal matrix elementcalculating unit configured to calculate a plurality of elements of anormal matrix for a least-squares method for calculating a displacementamount displaced from the preliminary alignment position, a pattern typeestimating unit configured to estimate a type of a pattern indicated bythe reference image, by using some of the plurality of elements of thenormal matrix, a displacement calculating unit configured to calculate asecond displacement amount displaced from the preliminary alignmentposition based on the least-squares method, by using a normal matrixobtained by deleting predetermined elements from the plurality ofelements of the normal matrix depending upon the type of the patternwhich has been estimated, a second SSD calculating unit configured tocalculate an SSD between the pixel value of the optical image and thepixel value of the reference image at a second position displaced fromthe preliminary alignment position by the second displacement amount, anSSD estimating unit configured to estimate which of the SSD at the firstposition and the SSD at the second position is smaller, a positioncorrecting unit configured to correct an alignment position between theoptical image and the reference image to a position where a smaller SSDas a result estimated by the SSD estimating unit is obtained, and acomparing unit configured to compare the optical image and the referenceimage whose alignment position has been corrected.

In accordance with another aspect of the present invention, a patterninspection apparatus includes an optical image acquiring unit configuredto acquire an optical image of pattern to be inspected on which apattern is formed, a reference image generating unit configured togenerate a reference image to be compared with the optical image, afirst SSD (Sum of Squared Difference) calculating unit configured tocalculate a first displacement amount from a preliminary alignmentposition between the optical image and the reference image to a firstposition where an SSD between a pixel value of the optical image and apixel value of the reference image is minimized, a normal matrix elementcalculating unit configured to calculate a plurality of elements of anormal matrix for a least-squares method for calculating a displacementamount displaced from the first position, a pattern type estimating unitconfigured to estimate a type of a pattern indicated by the referenceimage, by using some of the plurality of elements of the normal matrix,a displacement calculating unit configured to calculate a seconddisplacement amount displaced from the first position based on theleast-squares method, by using a normal matrix obtained by deletingpredetermined elements from the plurality of elements of the normalmatrix depending upon the type of the pattern which has been estimated,a second SSD calculating unit configured to calculate an SSD between thepixel value of the optical image and the pixel value of the referenceimage at a second position displaced from the first position by thesecond displacement amount, an SSD estimating unit configured toestimate which of the SSD at the first position and the SSD at thesecond position is smaller, a position correcting unit configured tocorrect an alignment position between the optical image and thereference image to a position where a smaller SSD as a result estimatedby the SSD estimating unit is obtained, and a comparing unit configuredto compare the optical image and the reference image whose alignmentposition has been corrected.

In accordance with another aspect of the present invention, an imagealignment method for aligning an optical image and a reference image foruse in a comparing inspection of pattern to be inspected on which apattern is formed, the method includes calculating a first displacementamount from a preliminary alignment position between the optical imageand the reference image to a first position where an SSD (Sum of SquaredDifference) between a pixel value of the optical image and a pixel valueof the reference image is minimized, calculating a plurality of elementsof a normal matrix for a least-squares method for calculating adisplacement amount displaced from the first position, estimating a typeof a pattern indicated by the reference image, by using some of theplurality of elements of the normal matrix, calculating a seconddisplacement amount displaced from the first position based on theleast-squares method, by using a normal matrix obtained by deletingpredetermined elements from the plurality of elements of the normalmatrix depending upon the type of the pattern which has been estimated,calculating an SSD between the pixel value of the optical image and thepixel value of the reference image at a second position displaced fromthe first position by the second displacement amount, estimating whichof the SSD at the first position and the SSD at the second position issmaller, and correcting an alignment position between the optical imageand the reference image to a position where a smaller SSD as a result ofthe estimating of SSD is obtained, to output a result of the correcting.

In accordance with another aspect of the present invention, adisplacement amount estimation method for estimating a displacementamount between an optical image and a reference image for use in acomparing inspection of pattern to be inspected on which a pattern isformed, by using a least-squares method, the displacement amountestimation method includes calculating a plurality of elements of anormal matrix for a least-squares method, based on a preliminaryalignment position between the optical image and the reference image,estimating a type of a pattern indicated by the reference image, byusing some of the plurality of elements of the normal matrix, andcalculating the displacement amount displaced from the preliminaryalignment position based on the least-squares method, by using a normalmatrix obtained by deleting predetermined elements from the plurality ofelements of the normal matrix depending upon the type of the patternwhich has been estimated, to output the displacement amount calculated.

In accordance with another aspect of the present invention, acomputer-readable recording medium with a program recorded causes acomputer to deploy the processes of storing an optical image and areference image for use in a comparing inspection of pattern to beinspected on which a pattern is formed, in a storage device, reading theoptical image and the reference image from the storage device, andcalculating a plurality of elements of a normal matrix for aleast-squares method, based on a preliminary alignment position betweenthe optical image and the reference image, estimating a type of apattern indicated by the reference image, by using some of the pluralityof elements of the normal matrix, and calculating a displacement amountdisplaced from the preliminary alignment position based on theleast-squares method, by using a normal matrix obtained by deletingpredetermined elements from the plurality of elements of the normalmatrix depending upon the type of the pattern which has been estimated,to output the displacement amount calculated.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic diagram showing the structure of patterninspection apparatus described in Embodiment 1;

FIG. 2 is a block diagram showing the structure of an alignment circuitdescribed in Embodiment 1;

FIG. 3 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 1;

FIG. 4 shows a drawing for explaining procedures of acquiring an opticalimage described in Embodiment 1;

FIGS. 5A, 5B, and 5C are shown for explaining dimension judgingdescribed in Embodiment 1;

FIG. 6 shows an example of a line & space pattern with an arbitraryangle described in Embodiment 1;

FIG. 7 is a diagram for explaining an angle range of the line & spacepattern with an arbitrary angle described in Embodiment 1;

FIGS. 8A and 8B are shown for explaining a method of deleting an elementof the line & space pattern with an arbitrary angle described inEmbodiment 1;

FIG. 9 shows an example of a two-dimensional pattern described inEmbodiment 2;

FIG. 10 is a block diagram showing the structure of an alignment circuitdescribed in Embodiment 2;

FIG. 11 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 2;

FIG. 12 shows a drawing for explaining weighting with two-dimensionallinear interpolation described in Embodiment 2;

FIG. 13 is a block diagram showing the structure of an alignment circuitdescribed in Embodiment 3;

FIG. 14 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 3;

FIG. 15 is a diagram for explaining an SSD calculation method describedin Embodiment 3;

FIG. 16 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 4; and

FIG. 17 is a diagram for explaining another method of acquiring anoptical image.

DETAILED DESCRIPTION OF THE INVENTION Embodiment 1

FIG. 1 is a schematic diagram showing the configuration of patterninspection apparatus described in Embodiment 1. In the figure, patterninspection apparatus 100 that inspects a defect of a substrate, such asa mask or a wafer on which a pattern is formed, serving as patternincludes an optical image acquiring unit 150 and a control systemcircuit 160. The optical image acquiring unit 150 includes an XYθ table102, a light source 103, a magnifying optical system 104, a photodiodearray 105, a sensor circuit 106, a laser length measurement system 122,an autoloader 130, and an illumination optical system 170. In thecontrol system circuit 160, a control calculator 110 serving as acomputer is connected, through a bus 120 serving as a data transmissionpath, to a position circuit 107, a comparing circuit 108, a referencecircuit 112 being an example of a reference image generating unit, analignment circuit 140, an autoloader control circuit 113, a tablecontrol circuit 114, a magnetic disk device 109, a magnetic tape device115, a flexible disk device (FD) 116, a CRT 117, a pattern monitor 118,and a printer 119. The XYθ table 102 is driven by an X-axis motor, aY-axis motor, and a θ-axis motor. In FIG. 1, only elements necessary forexplaining Embodiment 1 are described, and others are not described. Itshould be understood that other constituent elements generally necessaryfor the target workpiece inspection apparatus 100 are included.

FIG. 2 is a block diagram showing an example of the configuration of thealignment circuit in Embodiment 1. In the figure, the alignment circuit140 includes a reference data memory 302, a measurement data memory 304,a least-squares method displacement calculating circuit 322, acalculation data memory 330, and a position correcting circuit 350. Theleast-squares method displacement calculating circuit 322 includes anormal matrix element calculating unit 370, a dimension judging unit380, and a displacement calculating unit 390. The dimension judging unit380 includes a determinant absolute value calculating unit 382 and atrace calculating unit 384. The alignment circuit 140 receives referencedata from the reference circuit 112 and measurement data from theoptical image acquiring unit 150, performs the alignment of these itemsof data, and outputs the reference data and the measurement data to thecomparing circuit 108. The data etc. calculated in the alignment circuit140 is stored in the calculation data memory 330 as needed.

FIG. 3 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 1. In the figure, the target workpieceinspection method deploys a series of steps including an optical imageacquiring step (S102), a reference data generating step (S104), analignment step, and a comparing step (S404). As the alignment step beingone example of an image alignment method, a series of steps including astoring step (S202), a least-squares method displacement calculatingstep (S300), and a position correcting step (S402) is deployed.Moreover, in the least-squares method displacement calculating step(S300), a series of steps including a normal matrix element calculatingstep (S302), a dimension judging step (S304) and a calculating step(S306) of a displacement amount, an image transmission loss ratio (or“image strength fluctuation rate”) and a graylevel offset is executed.In FIG. 3, a solid line shows a flow of measurement data (opticalimage), and a dotted line shows a flow of reference data.

In S (step) 102, as the optical image acquiring step, the optical imageacquiring unit 150 acquires an optical image of a photomask 101 servingas pattern on which a figure indicated by figure data included in designdata is drawn based on the design data. More specifically, the opticalimage can be acquired as follows:

The photomask 101 serving as pattern to be inspected is placed on theXYθ table 102 which is movable in a horizontal direction and a rotatingdirection by the X-, Y-, and θ-axis motors. The pattern written on thephotomask 101 is irradiated with lights from the appropriate lightsource 103 arranged above the XYθ table 102. The photomask 101 servingas pattern is irradiated with a flux of light from the light source 103,through the illumination optical system 170. Below the photomask 101,the magnifying optical system 104, the photodiode array 105, and thesensor circuit 106 are arranged. The light transmitted through thephotomask 101 serving as pattern such as an exposure mask is focused onthe photodiode array 105 as an optical image, through the magnifyingoptical system 104 and enters the photodiode array 105.

FIG. 4 shows a diagram for explaining a procedure for acquiring anoptical image described in Embodiment 1. As shown in the figure, aregion to be inspected is virtually divided into a plurality ofstrip-like inspection stripes, each of which has a scanning swath W, inthe Y direction. To acquire an optical image, the movement of the XYθtable 102 is controlled so that each of the divided inspection stripescan be continuously scanned, while moving in the X direction. In thephotodiode array 105, images each having the scanning swath W as shownin FIG. 4 are continuously input. After images on the first inspectionstripe having been scanned, images each having the scanning swath W arecontinuously input while an image on the second inspection stripe ismoved in the reverse direction. When an image on the third inspectionstripe is to be acquired, the image is scanned while the image is movedin the direction reverse to the direction for scanning the image on thesecond inspection stripe, i.e., the same direction for scanning theimage on the first inspection stripe. Continuously acquiring images inthis manner makes it possible to reduce wasteful processing time.

The image of the pattern focused on the photodiode array 105 isphotoelectrically converted by the photodiode array 105. Furthermore,the electric image is A/D-converted (analog to digital converted) by thesensor circuit 106. In the photodiode array 105, a sensor such as a TDI(Time Delay Integration) sensor is arranged. The TDI sensor scans theimage of the pattern of the photomask 101 serving as pattern, bycontinuously moving the XYθ table 102 serving as a stage in the X-axisdirection. The light source 103, the magnifying optical system 104, thephotodiode array 105, and the sensor circuit 106 compose an inspectionoptical system having a large magnification.

The XYθ table 102 is driven by the table control circuit 114 under thecontrol of the control calculator 110. The XYθ table 102 can be moved bya drive system such as a three-axis (X-Y-θ) motor which drives the XYθtable 102 in the X direction, the Y direction, and the θ direction.

Measurement data (optical image) output from the sensor circuit 106 istransmitted to the alignment circuit 140 together with data which isoutput from the position circuit 107 and indicates the position of thephotomask 101 on the XYθ table 102. The measurement pattern data is, forexample, 8-bit unsigned data, and indicates a graylevel of brightness ofeach pixel. The measurement data is compared with each image data of512×512 pixels, for example.

Then, in step S104, as the reference data generating step, the referencecircuit 112 generates reference data (reference image) for comparingwith measurement data on the basis of design data of the photo mask 101serving as pattern to be inspected. The reference data to be compared isgenerated as image data of 512×512 pixels, for example, like themeasurement data.

The reference data herein is generated based on the design data in orderto execute a “die to database inspection”, but it does not restricted tothis. A “die to die inspection” can also be conducted, and in this case,reference data can be generated based on another measurement data(optical image) to be used for comparison.

Next, as the alignment step, aligning is performed for comparing themeasurement data and the reference data.

In step S202, as the storing step, reference data, for each 512×512pixels for example as stated above, is read using the control calculator110 and stored in the reference data memory 302. In the same manner,measurement data, for each 512×512 pixels for example, is read andstored in the measurement data memory 304. Next, a least-squares methoddisplacement calculation is performed. An amount of displacement neededfor alignment is herein calculated using the least-squares method beinga statistical method.

Supposing that a graylevel value (pixel value) of measurement dataserving as an optical image (actual image) is S (x, y), a displacementamount in the directions of X and Y of the graylevel value S(x, y) ofthe measurement data is (x₀, y₀), an image transmission loss ratio isε_(n) a graylevel offset is δ, and a pixel number is N+1, the equation(1) shown below can be obtained with respect to a graylevel value U(x,y) of reference data serving as a reference image.

$\begin{matrix}{{S( {x,y} )} = {{U( {{x - x_{0}},{y - y_{0}}} )} - {\sum\limits_{i = 0}^{N}\; {ɛ_{i}{U^{i + 1}( {{x - x_{0}},{y - y_{0}}} )}}} - \delta}} & (1)\end{matrix}$

Moreover, by linearization on supposition that the fluctuation amount issmall enough, the equation (2) shown below can be obtained.

$\begin{matrix}{{{U( {x,y} )} - {S( {x,y} )}} = {{x_{0}\frac{\partial U}{\partial x}} + {y_{0}\frac{\partial U}{\partial y}} + {\sum\limits_{i = 0}^{N}\; {ɛ_{i}U^{i + 1}}} + \delta}} & (2)\end{matrix}$

wherein ∂U/∂x is a partial differential (space differentiation) of U byx, and ∂U/∂y lay is a partial differential (space differentiation) of Uby y.

In step S302, as the normal matrix element calculating step, the normalmatrix element calculating unit 370 calculates a plurality of elementsof the normal matrix for the least-squares method, for calculating adisplacement amount (x₀, y₀) displaced from the preliminary alignmentposition between the measurement data and the reference data.Specifically, with respect to each pixel of a two-dimensional image, agraylevel value U (x, y) of reference data, a value (U-S) obtained bysubtracting a graylevel of the measurement data serving as an actualimage from a graylevel of the reference data, a value (∂U/∂x) obtainedby space differentiating the graylevel of the reference data in the Xdirection, and a value (∂U/∂y) obtained by space differentiating thegraylevel of the reference data in the Y direction are calculated toobtain each element of the equation (3) of the correlation matrix shownbelow. Moreover, as the preliminary alignment position, a positiontentatively in accordance in the data coordinate system can be used.

$\begin{matrix}{{\begin{pmatrix}{\sum\; ( \frac{\partial U}{\partial x} )^{2}} & {\sum\; {\frac{\partial U}{\partial x} \cdot \frac{\partial U}{\partial y}}} & {{\sum\; {\frac{\partial U}{\partial x} \cdot U}}\mspace{20mu}} & \cdots & {\sum\; {\frac{\partial U}{\partial x} \cdot U^{N}}} & {\sum\; \frac{\partial U}{\partial x}} \\{\sum\; {\frac{\partial U}{\partial y} \cdot \frac{\partial U}{\partial x}}} & {\sum\; ( \frac{\partial U}{\partial y} )^{2}} & {{\sum\; {\frac{\partial U}{\partial y} \cdot U}}\;} & \cdots & {\sum\; {\frac{\partial U}{\partial y} \cdot U^{N}}} & {\sum\; \frac{\partial U}{\partial y}} \\{\sum\; {U \cdot \frac{\partial U}{\partial x}}} & {\sum\; {U \cdot \frac{\partial U}{\partial y}}} & {\sum\; U^{2}} & \cdots & {\sum\; U^{N + 1}} & {\sum\; U} \\\vdots & \vdots & \vdots & \; & \vdots & \vdots \\{\sum{U^{N} \cdot \frac{\partial U}{\partial x}}} & {\sum{U^{N} \cdot \frac{\partial U}{\partial y}}} & {\sum\; U^{N + 1}} & \cdots & {\sum\; U^{2N}} & {\sum\; U^{N}} \\{\sum\; \frac{\partial U}{\partial x}} & {\sum\; \frac{\partial U}{\partial y}} & {{\sum\; U}\;} & \cdots & {\sum\; U^{N}} & 1\end{pmatrix}\begin{pmatrix}x_{0} \\y_{0} \\ɛ_{0} \\\vdots \\ɛ_{N} \\\delta\end{pmatrix}}\begin{pmatrix}{\frac{\partial U}{\partial x}( {U - S} )} \\{\frac{\partial U}{\partial y}( {U - S} )} \\{\sum\; {U( {U - S} )}} \\\vdots \\{\sum\; {U^{N}( {U - S} )}} \\{\sum\; ( {U - S} )}\end{pmatrix}} & (3)\end{matrix}$

By solving the simultaneous equations (3), the displacement amount (x₀,y₀) in the directions of X and Y, the image transmission loss ratios ε₀to ε_(N), and the graylevel off set δ can be obtained. In the patternwhose whole image range serves as a pattern (zero-dimensional pattern),such as a so-called no-structure pattern, the displacement amount (x₀,y₀) in the directions X and Y should intrinsically be indeterminate.However, there is actually a case that a wrong value has been calculatedas a solution, by influence of a noise etc. Moreover, for example, inthe pattern which spreads in one direction with an arbitrary angle(one-dimensional pattern), such as a so-called line and space pattern(line & space pattern), the displacement amount (x₀ or y₀) in thedirection X or Y should intrinsically be indeterminate. Actually, thereis a case that the solution becomes unstable.

Then, according to the present Embodiment, the type of a pattern to becompared is estimated. That is, it is estimated whether it is azero-dimensional pattern, such as an no-structure pattern, aone-dimensional pattern with an arbitrary angle, such as a line & spacepattern, or a two-dimensional pattern, such as a hole pattern and anL-shaped pattern. Then, solving is performed based on the methodcorresponding to the estimating result, and thereby an unstable solutioncan be eliminated.

In S304, as the dimension estimating step serving as an example of apattern type estimating unit, the dimension estimating unit 380estimates the type (dimension) of a pattern indicated by reference data,using some of a plurality of elements of the normal matrix shown in theequation (3). Specifically, the dimension is estimated by using thematrix (4) with two rows and two columns which is some of a plurality ofelements of the normal matrix shown in the equation (3).

$\begin{matrix}\begin{pmatrix}{\sum\; ( \frac{\partial U}{\partial x} )^{2}} & {\sum\; {\frac{\partial U}{\partial x}\frac{\partial U}{\partial y}}} \\{\sum\; {\frac{\partial U}{\partial x}\frac{\partial U}{\partial y}}} & {\sum\; ( \frac{\partial U}{\partial y} )^{2}}\end{pmatrix} & (4)\end{matrix}$

As shown in the matrix (4), the following are used as elements of the2×2 matrix: the total sum of differentiation values (∂U/∂x) (firstdifferentiation value) obtained by space differentiating a graylevelvalue of reference data in the X direction with respect to each pixel,the total sum of differentiation values (∂U/∂y) (second differentiationvalue) obtained by space differentiating a graylevel value of referencedata in the Y direction with respect to each pixel, the total sum ofsquared differentiation values (∂U/∂x), and the total sum of squareddifferentiation values (∂U/∂y). Then, the trace calculating unit 384calculates a trace of the 2×2 matrix shown in the matrix (4). Thedeterminant absolute value calculating unit 382 calculates an absolutevalue of the determinant of the 2×2 matrix shown in the matrix (4).

FIGS. 5A to 5C are shown for explaining dimension estimating describedin Embodiment 1. The estimating is performed using predeterminedthreshold values θ₁ and θ₂. For example, when reference data is aso-called no-structure pattern (zero-dimension), the absolute value of adeterminant of a matrix ideally becomes a value 0 (zero). When referencedata is a so-called line pattern (one-dimension), the trace ideallybecomes a value 0 (zero). In this case, the ideal value “0” (zero) canbe obtained when the reference data is generated from design data.However, there may be a case that when the reference data is generatedfrom measurement data, like the case of “die to die inspection”, theideal value “0” cannot be obtained because it is affected by a noiseetc. Then, in Embodiment 1, the estimating is performed by usingpredetermined threshold values θ₁ and θ₂ which are in consideration ofthe case that it is affected by the noise etc.

As shown in FIG. 5A, when the trace is smaller than the threshold valueθ₁ (first threshold value) and the absolute value of the matrixdeterminant is smaller than the threshold value θ₂ (second thresholdvalue), it is estimated that the type of a pattern is an no-structurepattern (zero-dimension) meaning the whole image being a pattern. Asshown in FIG. 5B, when the trace is greater than or equal to thethreshold value θ₁ and the absolute value of the matrix determinant issmaller than the threshold value θ₂, it is estimated that the type of apattern is a line & space pattern (one-dimension) with an arbitraryangle. As shown in FIG. 5C, when the trace is greater than or equal tothe threshold value θ₁ and the absolute value of the matrix determinantis greater than or equal to the threshold value θ₂, it is estimated tobe other pattern (two-dimension).

When estimated to be zero-dimension, since the differentiation value(∂U/∂x) and the differentiation value (∂U/∂y) intrinsically become “0”,the equation (3) of the correlation matrix can be rank deficient (ranknumber 1) as shown in the equation (5) by deleting the normal matrixelements including the above differentiation values.

$\begin{matrix}{{\begin{pmatrix}{\sum\; U^{2}} & \cdots & {\sum\; U^{N + 1}} & {\sum\; U} \\\vdots & \vdots & \vdots & \vdots \\{\sum\; U^{N + 1}} & \cdots & {\sum\; U^{2N}} & {\sum\; U^{N}} \\{\sum\; U} & \cdots & {\sum\; U^{N}} & 1\end{pmatrix}\begin{pmatrix}ɛ_{0} \\\vdots \\ɛ_{N} \\\delta\end{pmatrix}} = \begin{pmatrix}{\sum\; {U( {U - S} )}} \\\vdots \\{\sum\; {U^{N}( {U - S} )}} \\{\sum\; ( {U - S} )}\end{pmatrix}} & (5)\end{matrix}$

By performing the above, the elements related with the displacementamount (x₀, y₀) in the directions X and Y whose solution intrinsicallybecomes indeterminate are deleted from the simultaneous equations.Consequently, it becomes possible to prevent calculating and obtaining awrong value with respect to the displacement amount (x₀, y₀). The normalmatrix elements including the differentiation value (∂u/∂x) and thedifferentiation value (∂U/∂y) are deleted herein, but the same resultcan be obtained by defining the elements to be a value 0 (zero) insteadof deleting.

When estimated to be one-dimension, the normal matrix element includingselected one of the differentiation value (∂U/∂x) and thedifferentiation value (∂U/∂y) is deleted as mentioned below. That is,the equation (3) of the correlation matrix becomes rank deficient asshown in the equation (6) or the equation (7) (rank number 2).

$\begin{matrix}{{\begin{pmatrix}{\sum\; ( \frac{\partial U}{\partial x} )^{2}} & {\sum\; {\frac{\partial U}{\partial x}U}} & \cdots & {\sum\; {\frac{\partial U}{\partial x}U^{N}}} & {\sum\; \frac{\partial U}{\partial x}} \\{\sum\; {U\frac{\partial U}{\partial x}}} & {\sum\; U^{2}} & \cdots & {\sum\; U^{N + 1}} & {\sum\; U} \\\vdots & \vdots & \; & \vdots & \vdots \\{\sum\; {U^{N}\frac{\partial U}{\partial x}}} & {\sum\; U^{N + 1}} & \cdots & {\sum\; U^{2N}} & {\sum\; U^{N}} \\{\sum\; \frac{\partial U}{\partial x}} & {\sum\; U} & \cdots & {\sum\; U^{N}} & 1\end{pmatrix}\begin{pmatrix}x_{0} \\ɛ_{0} \\\vdots \\ɛ_{N} \\\delta\end{pmatrix}} = \begin{pmatrix}{\sum\; {\frac{\partial U}{\partial x}( {U - S} )}} \\{\sum\; {U( {U - S} )}} \\\vdots \\{\sum\; {U^{N}( {U - S} )}} \\{\sum\; ( {U - S} )}\end{pmatrix}} & (6)\end{matrix}$

Since the normal matrix elements including the differentiation value(∂U/∂y) are deleted in the equation (6), the equation (3) of thecorrelation matrix becomes rank deficient (rank number 2). As theresult, y₀ can be regarded to be indeterminate.

$\begin{matrix}{{\begin{pmatrix}{\sum\; ( \frac{\partial U}{\partial y} )^{2}} & {\sum\; {\frac{\partial U}{\partial y}U}} & \cdots & {\sum\; {\frac{\partial U}{\partial y}U^{N}}} & {\sum\; \frac{\partial U}{\partial y}} \\{\sum\; {U\frac{\partial U}{\partial y}}} & {\sum\; U^{2}} & \cdots & {\sum\; U^{N + 1}} & {\sum\; U} \\\vdots & \vdots & \; & \vdots & \vdots \\{\sum\; {U^{N}\frac{\partial U}{\partial y}}} & {\sum\; U^{N + 1}} & \cdots & {\sum\; U^{2N}} & {\sum\; U^{N}} \\{\sum\; \frac{\partial U}{\partial y}} & {\sum\; U} & \cdots & {\sum\; U^{N}} & 1\end{pmatrix}\begin{pmatrix}y_{0} \\ɛ_{0} \\\vdots \\ɛ_{N} \\\delta\end{pmatrix}} = \begin{pmatrix}{\sum\; {\frac{\partial U}{\partial y}( {U - S} )}} \\{\sum\; {U( {U - S} )}} \\\vdots \\{\sum\; {U^{N}( {U - S} )}} \\{\sum\; ( {U - S} )}\end{pmatrix}} & (7)\end{matrix}$

Since the normal matrix elements including the differentiation value(∂U/∂x) are deleted in the equation (7), the equation (3) of thecorrelation matrix becomes rank deficient (rank number 2). As theresult, x₀ can be regarded to be indeterminate.

FIG. 6 shows an example of the line & space pattern with an arbitraryangle described in Embodiment 1. In the figure, reference data 10composed of a line pattern 12 and a space pattern 14 is shown as anexample of the line & space pattern. The direction of the line & spacepattern does not necessarily correspond with the X-axis or the Y-axis ofa data coordinate system. The case of the line & space pattern (φ=90degrees) in the direction of Y will be explained as an example. First,the total sum of (∂U/∂x)² and the total sum of (∂U/∂y)², which are thenormal matrix elements of the equation (3), are compared. Then, sincethe total sum of (∂U/∂y)² can be ignored when compared with the totalsum of (∂U/∂x)², the normal matrix degenerates (rank deficient) bydeleting the term including (∂U/∂y) from the correlation equation (3) ofthe normal matrix. As a result, y₀ can be regarded to be indeterminate.Similarly, the case of the line & space pattern (φ=0 degree) in thedirection of X is now explained as an example. The total sum of (∂U/∂x)²and the total sum of (∂U/∂y)² which are the normal matrix elements ofthe equation (3) are compared. Then, since the total sum of (∂U/∂x)² canbe ignored when compared with the total sum of (∂U/∂y)², the normalmatrix degenerates (rank deficient) by deleting the term including(∂U/∂x) from the correlation equation (3) of the normal matrix. As aresult, x₀ can be regarded to be indeterminate. However, in the case ofthe line & space pattern with an arbitrary angle, a difference betweenthe total sum of (∂U/∂x)² and the total sum of (∂U/∂y)² is sometimes notlarge enough to be ignored. Therefore, it is difficult to estimate aterm to be deleted based on this method. Accordingly, in Embodiment 1,when estimated to be one-dimension as stated above, the element whichincludes either (∂U/∂x) or (∂U/∂y) is deleted by the following method.

FIG. 7 is a diagram for explaining an angle range of the line & spacepattern with an arbitrary angle described in Embodiment 1. In thefigure, the X direction is 0 degree and the Y direction is 90 degrees,for example. In the reference data in FIG. 7, the case of the angle φ ofthe line & space pattern with an arbitrary angle being greater than orequal to 0 degree and less than 45 degrees or the case of the angle φbeing greater than or equal to 135 degrees and less than or equal to 180degrees is defined as a range A. The case of the angle φ of the line &space pattern with an arbitrary angle being greater than or equal to 45degrees and less than 135 degrees is defined as a range B. In eachrange, the total sum of (∂U/∂x) 2 and the total sum of (∂U/∂y)² whichare elements of the normal matrix of the equation (3) are compared.

FIGS. 8A and 8B are shown for explaining a method of deleting an elementof the line & space pattern with an arbitrary angle described inEmbodiment 1. When comparing the total sum of (∂U/∂x)² and the total sumof (∂U/∂y)² which are diagonal components, the total sum of (∂U/∂x)² issmaller than the total sum of (∂U/∂y)² in the range A, as shown in FIG.8A. Then, in the case of a line & space pattern with an arbitrary anglecorresponding to the range A in the reference data, the term including(∂U/∂x) is deleted from the correlation equation (3) of the normalmatrix. This makes the normal matrix degenerate (rank deficient). As aresult, x₀ can be regarded to be indeterminate. On the other hand, inthe range B of the reference data, the total sum of (∂U/∂y)² is smallerthan the total sum of (∂U/∂x)² as shown in FIG. 8B. Then, in the case ofa line & space pattern with an arbitrary angle corresponding to therange B, the term including (∂U/∂y) is deleted from the correlationequation (3) of the normal matrix. This makes the normal matrixdegenerate (rank deficient). As a result, y₀ can be regarded to beindeterminate.

That is, even in the case of the line & space pattern with an arbitraryangle, the above configuration makes it possible to delete an elementincluding either one of (∂U/∂x) and (∂U/∂y) which are preferably to beintrinsically indeterminate.

When estimated to be two-dimension, since a solution being indeterminateis not generated, it is just necessary to solve the simultaneousequations of the correlation matrix shown in the equation (3) withoutdeleting any element of the normal matrix (rank number 3).

As mentioned above, a correlation matrix corresponding to the patterntype can be obtained by estimating a dimension by using a trace and anabsolute value of a determinant of matrix.

In step S306, as the calculating step of a displacement amount, an imagetransmission loss ratio, and a graylevel offset, which serves an exampleof the displacement amount calculating step, the displacementcalculating unit 390 solves simultaneous equations of the normal matrixshown in the equation (5), (6), or (7) in which a predetermined elementis deleted from a plurality of elements of the normal matrix, dependingupon a pattern type respectively estimated, or solves the normal matrixshown in the equation (3) in which no element is deleted. By performingsuch a simultaneous equation solution step, a displacement amount (x₀,y₀) displaced from a preliminary alignment position obtained by theleast-squares method, an image transmission loss ratio ε₀ to ε_(N) and agraylevel offset δ can be calculated depending upon the type of apattern. In the above explanation, description relating to each value ofthe matrix is partially omitted in the equations (3), (5), (6) and (7).

The above configuration makes it possible to achieve a solutiondepending upon a solution method suitable for each pattern type.Therefore, it becomes possible to obtain a displacement amount (x₀, y₀),an image transmission loss ratio ε₀ to ε_(N) and a graylevel offset δfor alignment to the best position based on the least-squares method.Then, in calculating these solutions, it is possible to eliminate asolution which should intrinsically be indeterminate.

In step S402, as the position correcting step, the position correctingcircuit 350 (an example of a position correcting unit) corrects analignment position between measurement data and reference data from thepreliminary alignment position to the position obtained by performingdisplacing by the displacement amount from the preliminary alignmentposition. In the case x₀ or y₀ being an indeterminate solution, thepreliminary alignment position can be used without changing, withrespect to the indeterminate solution. It is also preferable for theposition correcting circuit 350 to correct the image strength of eachpixel of reference data using the image transmission loss ratio ε₀ toε_(N) calculated by the least-squares method displacement calculatingcircuit 322. The result of the correcting is output to the comparingcircuit 108.

In step S404, as the comparing step, the comparing circuit 108 aligns,by means of the position alignment circuit 140, the measurement dataserving as a pattern image to be inspected generated by the sensorcircuit 106 on the basis of the transfer image obtained from the photomask 101, and the reference data serving as an inspection standardpattern image generated by the reference circuit 112, and takes in boththe data. Then, the comparing circuit 108 compares them, namely thetaken measurement data and reference data, with each other according toa predetermined algorithm, and estimates whether there is a defect ornot. The comparing circuit 108 outputs the result of the comparing.Thus, by performing a data comparison through such a highly precisealignment, it is possible to prevent a false detection of a defect andto decrease pseudo defects, thereby performing a highly preciseinspection.

Embodiment 2

A method of correcting a local displacement in a frame will be explainedin Embodiment 2.

FIG. 9 shows an example of a two-dimensional pattern described inEmbodiment 2. For example, the case where a local displacement occurs ina part of the two-dimensional pattern shown in FIG. 9 will be explained.In the least-squares method, it is preferable to divide such a frame bydot lines as shown in FIG. 9 for example, to calculate a displacementamount and an image transmission loss ratio for each of the dividedareas, and to estimate the displacement amount and the imagetransmission loss ratio, respectively. In such a case, when the dividedareas are configured by a zero dimensional pattern, normal matrices canbe degenerated to make x₀ and y₀ become indeterminate, as explained inEmbodiment 1. Similarly, when the divided areas are configured by aone-dimensional pattern, regular matrices can be degeneratedrespectively to make either x₀ or y₀ become indeterminate, as explainedin Embodiment 1. Consequently, it is possible to perform a highlyprecise displacement correction by composing the respective dividedareas.

Further, in the least-squares method, when dividing the above-mentionedframe by dot lines as shown in FIG. 8, it is also preferable to dividethe frame by weighting, to calculate a displacement amount for each ofthe divided areas, and to estimate the displacement amount and the imagetransmission loss ratio, respectively.

FIG. 10 is a block diagram showing the structure of the alignmentcircuit in Embodiment 2. In the figure, the alignment circuit 140further includes a weighting unit 360 in the least-squares methoddisplacement calculating circuit 322 in addition to the structure shownin FIG. 2. Namely, FIG. 10 is the same as FIG. 2 except for theweighting unit 360 added.

FIG. 11 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 2. In the figure, a weighting step (S301) isadded before the normal matrix element calculating step (S302) in theleast-squares method displacement calculating step (S300) in addition tothe steps shown in FIG. 3. Namely, FIG. 11 is the same as FIG. 3 exceptfor the weighting step S301 added. Furthermore, each step of theapparatus structure and the target workpiece inspection method or theimage alignment method in Embodiment 2 is the same as that of Embodiment1.

In step S301, as the weighting step, the weighting unit 360 calculates aweighting factor for performing weighting by one-dimensional linearinterpolation. The weighting unit 360 calculates a value multiplied bythe weighting factor with respect to each graylevel value of referencedata.

FIG. 12 is a drawing for explaining weighting by one-dimensional linearinterpolation in Embodiment 2. The case where weighting is performed atneighboring four points will now be explained. For example, with regardto certain image data P, when a pixel P is interpolated by use of imagedata of four points (0, 1, 2, 3) therearound, the following equation (8)can be defined.

$\begin{matrix}( \begin{matrix}{U_{(0)} = {{( {1 - x} )( {1 - y} )U} = {W_{(0)}U}}} \\{U_{(1)} = {{{x( {1 - y} )}U} = {W_{(1)}U}}} \\{U_{(2)} = {{xyU} = {W_{(2)}U}}} \\{U_{(3)} = {{( {1 - x} ){yU}} = {W_{(3)}U}}}\end{matrix}  & (8)\end{matrix}$

That is, as shown in a method (8), a graylevel value U₍₀₎ of referencedata can be expressed by U(0)=(1−x) (1−y) U. When (1−x) (1−y) expressesa weighting factor W₍₀₎, the graylevel value can be expressed asU₍₀₎=W₍₀₎ U. Moreover, reference data U₍₁₎ can be expressed byU₍₁₎=x(1−y) U. When x(1−y) expresses a weighting factor W₍₁₎, thereference data can be expressed as U₍₁₎=W₍₁₎ U. Moreover, reference dataU₍₂₎ can be expressed by U₍₂₎=xyU. When xy expresses a weighting factorW₍₂₎, the reference data can be expressed as U₍₂₎=W₍₂₎ U. Moreover,reference data U₍₃₎ can be expressed by U₍₃₎=(1−x) yU. When (1−x)yexpresses a weighting factor W₍₃₎, the reference data can be expressedas U₍₃₎=W₍₃₎ U. Thus, if a correlation matrix is calculated by use ofthe weighted reference data U₍₀₎, reference data U₍₁₎, reference dataU₍₂₎, and reference data U₍₃₎, the correlation matrix can be shown asthe equation (9).

$\begin{matrix}{{\begin{pmatrix}a & b & c & d \\e & f & g & h \\i & j & k & l \\m & n & p & q\end{pmatrix}\begin{pmatrix}A \\B \\C \\D\end{pmatrix}} = \begin{pmatrix}A^{\prime} \\B^{\prime} \\C^{\prime} \\D^{\prime}\end{pmatrix}} & (9)\end{matrix}$

where each element a to d in the equation (9) can be expressed as thefollowing equation (10).

$\quad\begin{matrix}\{ \begin{matrix}{{a = \begin{pmatrix}{\sum\; ( \frac{\partial U_{(0)}}{\partial x} )^{2}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial x}\frac{\partial U_{(3)}}{\partial x}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial x}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; ( \frac{\partial U_{(3)}}{\partial x} )^{2}}\end{pmatrix}},} & {c = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial x}U_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial x}U_{(3)}^{N}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial x}U_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial x}U_{(3)}^{N}}}\end{pmatrix}} \\{{b = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial x}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial x}\frac{\partial U_{(3)}}{\partial y}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial x}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial x}\frac{\partial U_{(3)}}{\partial y}}}\end{pmatrix}},} & {d = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial x}W_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial x}W_{(3)}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial x}W_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial x}W_{(3)}}}\end{pmatrix}}\end{matrix}  & (10)\end{matrix}$

Moreover, each element e to h in the equation (9) can be expressed asthe following equation (11).

$\quad\begin{matrix}\{ \begin{matrix}{{e = \begin{pmatrix}{\sum{\frac{\partial U_{(0)}}{\partial y}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial y}\frac{\partial U_{(3)}}{\partial x}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial y}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial y}\frac{\partial U_{(3)}}{\partial x}}}\end{pmatrix}},} & {g = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial y}U_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial y}U_{(3)}^{N}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial y}U_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial y}U_{(3)}^{N}}}\end{pmatrix}} \\{{f = \begin{pmatrix}{\sum\; ( \frac{\partial U_{(0)}}{\partial y} )^{2}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial y}\frac{\partial U_{(3)}}{\partial y}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial y}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; ( \frac{\partial U_{(3)}}{\partial y} )^{2}}\end{pmatrix}},} & {h = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial y}W_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(0)}}{\partial y}W_{(3)}}} \\\vdots & ⋰ & \vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial y}W_{(0)}}} & \cdots & {\sum\; {\frac{\partial U_{(3)}}{\partial y}W_{(3)}}}\end{pmatrix}}\end{matrix}  & (11)\end{matrix}$

Moreover, each element i to l in the equation (9) can be expressed asthe following equation (12).

$\quad\begin{matrix}\{ \begin{matrix}{{i = \begin{pmatrix}{\sum\; {U_{(0)}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {U_{(0)}\frac{\partial U_{(3)}}{\partial x}}} \\\vdots & ⋰ & \vdots \\{\sum\; {U_{(3)}^{N}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {U_{(3)}^{N}\frac{\partial U_{(3)}}{\partial x}}}\end{pmatrix}},} & {k = \begin{pmatrix}{\sum\; U_{(0)}^{2}} & \cdots & {\sum\; {U_{(0)}U_{(3)}^{N}}} \\\vdots & ⋰ & \vdots \\{\sum\; {U_{(3)}^{N}U_{(0)}}} & \cdots & {\sum\; U_{(3)}^{2N}}\end{pmatrix}} \\{{j = \begin{pmatrix}{\sum\; {U_{(0)}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {U_{(0)}\frac{\partial U_{(3)}}{\partial y}}} \\\vdots & ⋰ & \vdots \\{\sum\; {U_{(3)}^{N}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {U_{(3)}^{N}\frac{\partial U_{(3)}}{\partial y}}}\end{pmatrix}},} & {l = \begin{pmatrix}{\sum\; {U_{(0)}W_{(0)}}} & \cdots & {\sum\; {U_{(0)}W_{(3)}}} \\\vdots & ⋰ & \vdots \\{\sum\; {U_{(3)}W_{(0)}}} & \cdots & {\sum\; {U_{(3)}W_{(3)}}}\end{pmatrix}}\end{matrix}  & (12)\end{matrix}$

Moreover, each element m to q in the equation (9) can be expressed asthe following equation (13).

$\begin{matrix}{\quad\{ \begin{matrix}{{m = \begin{pmatrix}{\sum\; {W_{(0)}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {W_{(0)}\frac{\partial U_{(3)}}{\partial x}}} \\\vdots & ⋰ & \vdots \\{\sum{W_{(3)}^{N}\frac{\partial U_{(0)}}{\partial x}}} & \cdots & {\sum\; {W_{(3)}^{N}\frac{\partial U_{(3)}}{\partial x}}}\end{pmatrix}},} & {p = \begin{pmatrix}{\sum\; {W_{(0)}U_{(0)}}} & \cdots & {\sum\; {W_{(0)}U_{(3)}^{N}}} \\\vdots & ⋰ & \vdots \\{\sum\; {W_{(3)}^{N}U_{(0)}}} & \cdots & {\sum\; {W_{(3)}U_{(3)}}}\end{pmatrix}} \\{{n = \begin{pmatrix}{\sum\; {W_{(0)}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {W_{(0)}\frac{\partial U_{(3)}}{\partial y}}} \\\vdots & ⋰ & \vdots \\{\sum\; {W_{(3)}\frac{\partial U_{(0)}}{\partial y}}} & \cdots & {\sum\; {W_{(3)}\frac{\partial U_{(3)}}{\partial y}}}\end{pmatrix}},} & {q = \begin{pmatrix}{\sum\; W_{(0)}^{2}} & \cdots & {\sum\; {W_{(0)}W_{(3)}}} \\\vdots & ⋰ & \vdots \\{\sum\; {W_{(3)}W_{(0)}}} & \cdots & {\sum\; W_{(3)}^{2}}\end{pmatrix}}\end{matrix} } & (13)\end{matrix}$

Moreover, each element A to D in the equation (9) can be expressed asthe following equation (14).

$\begin{matrix}\{ \begin{matrix}{{A = \begin{pmatrix}x_{0_{(0)}} \\\vdots \\x_{0_{(3)}}\end{pmatrix}},} & {C = \begin{pmatrix}ɛ_{0_{(0)}} \\\vdots \\ɛ_{N_{(3)}}\end{pmatrix}} \\{{B = \begin{pmatrix}y_{0_{(0)}} \\\vdots \\y_{0_{(3)}}\end{pmatrix}},} & {D = \begin{pmatrix}\delta_{(0)} \\\vdots \\\delta_{(3)}\end{pmatrix}}\end{matrix}  & (14)\end{matrix}$

Moreover, each element A′ to D′ in the equation (9) can be expressed asthe following equation (15).

$\begin{matrix}{\quad\{ \begin{matrix}{{A^{\prime} = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial x}( {U_{(0)} - S_{(0)}} )}} \\\vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial x}( {U_{(3)} - S_{(3)}} )}}\end{pmatrix}},} & {C^{\prime} = \begin{pmatrix}{\sum\; {U_{(0)}( {U_{(0)} - S_{(0)}} )}} \\\vdots \\{\sum\; {U_{(3)}^{N}( {U_{(3)} - S_{(3)}} )}}\end{pmatrix}} \\{{B^{\prime} = \begin{pmatrix}{\sum\; {\frac{\partial U_{(0)}}{\partial y}( {U_{(0)} - S_{(0)}} )}} \\\vdots \\{\sum\; {\frac{\partial U_{(3)}}{\partial y}( {U_{(3)} - S_{(3)}} )}}\end{pmatrix}},} & {D^{\prime} = \begin{pmatrix}{\sum\; {W_{(0)}( {U_{(0)} - S_{(0)}} )}} \\\vdots \\{\sum\; {W_{(3)}( {U_{(3)} - S_{(3)}} )}}\end{pmatrix}}\end{matrix} } & (15)\end{matrix}$

As mentioned above, when weighting is performed at neighboring fourpoints, the correlation matrix equations shown in equations (9) to (15)are solved to obtain 4(N+4) variables, such as image transmission lossratios ε₀₍₀₎ to ε_(N(3)), displacement amounts x₀₍₀₎, x₀₍₁₎, x₀₍₂₎,x₀₍₃₎, y₀₍₀₎, y₀₍₁₎, y₀₍₂₎, and y₀₍₃₎, and graylevel offsets δ₍₀₎ toδ₍₃₎. When indeterminate solutions exist, it should be understood thatthe number of the indeterminate solutions are to be excluded. Positioncorrection can be performed by use of such values to correct a localdisplacement and the like in the frame. Causes of the local displacementmay include (1) meandering of the XY stage, (2) a pixel size differencebetween an actual image and a reference image, and (3) a pixel sizedifference between image scanning elements. The weighting method is notlimited to the one using neighboring 4 points. It is also preferable toadopt bicubic interpolation using 16 points or the like. In theequations (1) to (15), description of values of the matrix is partlyomitted.

Embodiment 3

In Embodiments 1 and 2 mentioned above, as the alignment method, adisplacement amount calculated by the least-squares method displacementcalculation in which dimension is estimated and a correlation matrix ismade to be rank deficient according to a pattern type is used. InEmbodiment 3, in addition to the least-squares method mentioned above,an alignment method where a SSD (Sum of the Squared Difference) methodis combined with the least-squares method will be explained.Specifically, the case where a sub-pixel unit SSD calculation and theleast-squares method displacement calculation are performed in parallelwill be described in Embodiment 3.

FIG. 13 is a block diagram showing an example of the structure of thealignment circuit described in Embodiment 3. In the figure, thealignment circuit 140 includes the reference data memory 302, themeasurement data memory 304, a pixel unit SSD calculating circuit 310, asub-pixel unit SSD calculating circuit 320, the least-squares methoddisplacement calculating circuit 322, the calculation data memory 330,an SSD (or “Residual Sum of Square (RSS)”) calculating circuit 324, aestimating circuit 340, and a position correcting circuit 350. Theposition alignment circuit 140 receives reference data from thereference circuit 112 and measurement data from the optical imageacquiring unit 150, performs the alignment of these items of data, andoutputs the reference data and the measurement data to the comparingcircuit 108. Each of the configuration in the least-squares methoddisplacement calculating circuit 322 is the same as that of FIG. 2 orFIG. 10.

FIG. 14 is a flowchart showing main steps of pattern inspection methoddescribed in Embodiment 3. In the figure, the target workpieceinspection method executes a series of steps including the optical imageacquiring step (S102), the reference data generating step (S104), thealignment step, and the comparing step (S404). As the alignment stepbeing one example of the image alignment method, a series of stepsincluding the storing step (S202), a pixel unit SSD calculating step(S204), a sub-pixel unit SSD calculating step (S206), the least-squaresmethod displacement calculating step (S300), a SSD calculating step(S308), a estimating step (S400), and the position correcting step(S402) are executed. In FIG. 14 as well as FIG. 3, solid lines show theflow of measurement data (optical image), and dotted lines show the flowof reference data.

Each of the steps from S102 to S202 is the same as that in Embodiment 1or Embodiment 2.

In step S204, as the pixel unit SSD calculating step, the pixel unit SSDcalculating circuit 310 serving as one example of an SSD calculatingunit calculates a displacement amount from a first preliminary alignmentposition to a position where the SSD between a pixel value of themeasurement data and a pixel value of the reference data becomes theminimum, by performing shifting in a pixel unit from the firstpreliminary alignment position. A position which should be tentativelyin accordance in the data coordinate system can be used as the firstpreliminary alignment position, which is the same as Embodiment 1 or 2.

FIG. 15 is a diagram for explaining an SSD calculation method describedin Embodiment 3. First, the pixel unit SSD calculating circuit 310 readsreference data of an image area of a predetermined size (frame) servingas the unit of comparing process, from the reference data memory 302 onthe basis of positional information from the position circuit 107. Atthis moment, the pixel unit SSD calculating circuit 310 respectivelygenerates images (displaced images), shifted in parallel in units ofpixels, with respect to the reference data of such a frame. In FIG. 15,the generated images are shown as data 1, data 2, . . . data n. Themeasurement data and the reference data in the frame are compared. Forexample, it is preferable to make an area of 512×512 pixels as oneframe. Between each item of plural reference data displaced in units ofpixels and the measurement data of the area of the same size read fromthe measurement data memory 304, the SSD is calculated. The SSD isobtained by summing the squared residual between each pixel value of thereference data and each pixel value of the measurement data. Then, theSSD of each of the plural reference data is calculated, and the minimumvalue of the SSD is calculated. The measurement data and the referencedata are aligned to a position where the minimum value is obtained. Inthis manner, it is possible to make alignment to the position where themeasurement data and the reference data are positioned closest whenshifted in parallel in x and y directions in units of pixels. Such aposition is made as a second preliminary alignment position, and adetailed alignment is performed hereinafter.

In step S206, as the sub-pixel image unit SSD calculating step, thesub-pixel unit SSD calculating circuit 320 serving as one example of theSSD calculating unit performs displacing in units of pixels from thesecond preliminary alignment position between the measurement data andthe reference data, which is preliminarily aligned in the pixel unit SSDcalculating step, to the position (first position) where the SSD betweena pixel value of the measurement data and a pixel value of the referencedata is minimized, and calculates the displacement amount (firstdisplacement amount) from the second preliminary alignment position tothe first position.

The sub-pixel unit SSD calculating method is the same as the contentsexplained with reference to FIG. 15. On the basis of the secondpreliminary alignment position, images (displaced images) shifted inparallel in units of sub-pixels are respectively generated with respectto the reference data of the size of the area to be compared. In FIG.15, the generated images are shown as data 1, data 2, . . . data n. Forexample, as sub-pixels, ⅛, 1/16, 1/32 and the like of one pixel are madeinto units. For example, when ⅛ of one pixel is made as the unit of thesub-pixel, the reference data of areas of a predetermined size displacedby +⅛ pixel, + 2/8 pixel, ±⅜ pixel, + 4/8 pixel, +⅝ pixel, ± 6/8 pixel,and +⅞ pixel in the X direction and the Y direction, respectively, andthe reference data with the displacement amount of 0 are generated. Thatis, 256 kinds of reference data in a combination of 16 ways in the Xdirection and 16 ways in the Y direction are generated. Then, the SSD iscalculated between the respective reference data and the respectivemeasurement data. The SSD is obtained by summing the squared residualbetween each pixel value of the reference data and each pixel value ofthe measurement data. Then, the SSD of each of the plural reference datais calculated, and the minimum value of the SSD is calculated. In thismanner, it is possible to obtain the displacement amount to the positionwhere the minimum value is obtained. The data such as the displacementamount which has been set and the calculated SSD are stored in thecalculation data memory 330. In this way, it is possible to obtain thedisplacement amount (x₀, y₀) for alignment of the measurement data andthe reference data to the position where they are positioned closestwhen shifted in parallel in the X and Y directions in units ofsub-pixels.

In step S300, as the least-squares method displacement calculating step,the least-squares method displacement calculating circuit 322 serving asone example of a least-squares method calculating unit calculates adisplacement amount (second displacement amount) based on theleast-squares method displaced from the above-mentioned secondpreliminary alignment position between the measurement data and thereference data. The contents of the displacement calculation based onthe least-squares method is the same as those in Embodiment 1 orEmbodiment 2. That is, after performing a dimension judging and makingthe correlation matrix be rank deficient depending upon the type of apattern, the displacement amount etc. is obtained by a displacementamount calculation based on the least-squares method. In other words,according to the method of Embodiment 1, it is possible to obtain 4(N+4)variables, such as image transmission loss ratios ε₀ to ε_(N),displacement amounts x₀, and y₀, and a graylevel offset δ. According tothe method of Embodiment 2, it is possible to obtain 4(N+4) variables,such as image transmission loss ratios ε₀₍₀₎ to ε_(N(3)), displacementamounts x₀₍₀₎, x₀₍₁₎, x₀₍₂₎, x₀₍₃₎, y₀₍₀₎, y₀₍₁₎, y₀₍₂₎, and y₀₍₃₎, andgraylevel offsets δ₍₀₎ to δ₍₃₎. In the case where dividing is performedby weighting as the method of Embodiment 2, the displacement amount (x₀,y₀) for alignment to the position where the measurement data and thereference data are positioned closest can be obtained by composing. Inthe methods of Embodiments 1 and 2, if indeterminate solutions exist,the indeterminate solutions are excluded as a matter of course.

In step S308, as the SSD calculating step, the SSD calculating circuit324 serving as one example of the SSD calculating unit calculates theSSD between a pixel value of the measurement data and a pixel value ofthe reference data at the position (x-x₀, y-y₀) (second position)displaced by the displacement amount (x₀, y₀) calculated by theleast-squares method displacement calculating circuit 322 from theabove-mentioned preliminary alignment position between the measurementdata and the reference data.

In step S400, as the estimating step, the estimating circuit 340 servingas an example of a estimating unit estimates which of the SSD at thefirst position and the SSD at the second position is smaller. That is,the estimating circuit 340 estimates which of the minimum SSD obtainedas the result of the calculation by the sub-pixel unit SSD calculatingcircuit 320 and the SSD obtained as the result of the calculation by theSSD calculating circuit 324 is smaller.

In step S402, as the position correcting step, the position correctingcircuit 350 serving as one example of the position correcting unitcorrects the alignment position between the measurement data and thereference data to a position where the smaller SSD determined by theestimating circuit 340 is obtained. Further, it is also preferable thatthe position correcting circuit 350 corrects the image graylevel of eachpixel of the reference data by use of the image transmission loss ratioE calculated by the least-squares method displacement calculatingcircuit 322. For example, not only when the estimating circuit 340adopts the result calculated by the SSD calculating circuit 324, butalso when the estimating circuit 340 adopts the result calculated by thesub-pixel unit SSD calculating circuit 320, the image graylevel of eachpixel of the reference data is preferably corrected by use of the imagetransmission loss ratio ε calculated by the least-squares methoddisplacement calculating circuit 322.

Herein, the type of patterns suitable for the SSD method or theleast-squares method is different. For example, the SSD method is suitedfor aligning patterns of sparse figure density. On the other hand, theleast-squares method is suited for aligning patterns of dense figuredensity. For this reason, with the configuration as explained in thepresent Embodiment, the SSD of the least-squares method is compared withthe minimum SSD of the SSD in units of sub-pixels, and the correctingmethod with the smaller SSD between the SSD of the least-squares methodand the SSD in units of sub-pixels is adopted, so that better resultsare expected than those in a case of correction made singly by each ofthe methods.

More specifically, in the case of an image of a sparse pattern, thecalculation by the least-squares method may become unstable, and thus,alignment by the SSD is adopted in that case. A parallel use of the SSDmethod and the least-squares method makes it possible to stably correcteven such a sparse pattern.

By correcting the reference data serving as a reference image or themeasurement data serving as an optical image (actual image) by use ofsuch a value, it is possible to make the measurement data and thereference data further closer to each other. As a result, it is possibleto prevent a false detection in inspecting defects, and to increase thepractical sensitivity. As mentioned above, by simply correcting thedisplacement between the reference image and the actual image and theimage transmission loss, a highly sensitive inspection can be achieved.

Herein, the object to be compared in the estimating step is not limitedto the SSD, but the sum of the p-th power of a residual wherein p is apositive number may be adopted generally. The SSD corresponds to thecase of P=2. In other words, a position correction by the SSD isperformed in parallel with the least-squares method, and the sum of thep-th power (p is a positive number) of the residual absolute value ofthe actual image and the corrected reference image is calculated in therespective cases of the correction by the least-squares method and thecorrection by the SSD method, both of the values are compared with eachother, and a correction method in which the sum of the p-th power of theresidual absolute value becomes minimum may be selected. Then, theresult of the correcting is output to the comparing circuit 108.

In step S404, as the comparing step, the comparing circuit 108 aligns,by means of the position alignment circuit 140, the measurement dataserving as a pattern image to be inspected generated by the sensorcircuit 106 on the basis of the transfer image obtained from the photomask 101, and the reference data serving as an inspection standardpattern image generated by the reference circuit 112, and then takes inboth the data. Then, the comparing circuit 108 compares them, namely thetaken measurement data and reference data, with each other according toa predetermined algorithm, and estimates whether there is a defect ornot. The comparing circuit 108 outputs the result of the comparing.Thus, by performing a data comparison through such a highly precisealignment, it is possible to prevent a false detection of a defect andto decrease pseudo defects, thereby performing a highly preciseinspection.

Embodiment 4

In Embodiment 3, the case where a sub-pixel unit SSD calculation isperformed in parallel with a least-squares method displacementcalculation has been explained. In Embodiment 4, the case where thesub-pixel unit SSD calculation is performed in series with theleast-squares method displacement calculation will now be explained. Theapparatus structure described in Embodiment 4 is the same as that inEmbodiment 3. Moreover, each configuration in the least-squares methoddisplacement calculating circuit 322 is the same as that in Embodiment3, and the same as that of FIG. 2 or FIG. 10.

FIG. 16 is a flowchart showing main steps of the target workpieceinspection method described in Embodiment 4.

In FIG. 14, the sub-pixel unit SSD calculating step (S206) is performedin parallel with a combination of the least-squares method displacementcalculating step (S300) and the SSD calculating step (S308). In FIG. 16,after the sub-pixel unit SSD calculating step (S206), the least-squaresmethod displacement calculating step (S300) and the SSD calculating step(S308) are performed, and others are the same as those of FIG. 14. InFIG. 16, as well as FIG. 3, solid lines show the flow of measurementdata (optical image), and dotted lines show the flow of reference data.

Each of the steps from S104 to S206 in FIG. 16 is the same as that ofEmbodiment 3.

In step S300 in FIG. 16, as the least-squares method displacementcalculating step, the least-squares method displacement calculatingcircuit 322 serving as one example of the least-squares methodcalculating unit calculates a displacement amount (second displacementamount) based on the least-squares method displaced from the firstposition of the measurement data and the reference data. That is,further alignment is performed for the position (first position) wherethe minimum SSD obtained as a result of calculating by the sub-pixelunit SSD calculating circuit 320. Thus, by performing further alignmentfrom the position obtained by the sub-pixel unit SSD calculation, highlyprecise alignment can be performed. The contents of the displacementamount calculation based on the least-squares method are the same asthose in Embodiment 1 or 2, that is, the same as those in Embodiment 3.

In step S308 in FIG. 16, as the SSD calculating step, the SSDcalculating circuit 324 serving as one example of the SSD calculatingunit calculates an SSD between the pixel value of the measurement dataand the pixel value of the reference data at the position (x-x₀, y-y₀)(second position) obtained by displaced by the displacement amount (x₀,y₀) calculated by the least-squares method displacement calculatingcircuit 322 from the above-mentioned first position of the measurementdata and the reference data.

In step S400 of FIG. 16, as the estimating step, the estimating circuit340 serving as one example of the estimating unit estimates which of theSSD at the first position and the SSD at the second position is smaller.That is, the judgment circuit 340 estimates which of the minimum SSDobtained as the result of the calculation by the sub-pixel unit SSDcalculating circuit 320 and the SSD obtained as the result of thecalculation by the SSD calculating circuit 324 is smaller. The followingis the same as Embodiment 3.

As mentioned above, by performing alignment from the position obtainedby the sub-pixel unit SSD calculation, it is possible to further enhancethe precision of Embodiment 3.

In each of the Embodiments mentioned above, with the configuration ofFIG. 1, the photodiode array 105 which simultaneously inputs beamscorresponding to the number of pixels (for example, 2048 pixels) of thescanning swath W is employed, but not limited to it. FIG. 17 is adiagram for explaining another method for acquiring an optical image. Asshown in FIG. 17, an alternative method may be used in which, while theXYθ table 102 is transferred at a constant speed in the X direction, alaser scanning optical device (not shown) scans with a laser beam in theY direction at every time when movement of a predetermined pitch isdetected by a laser interferometer, and transmitted light is detected toacquire a two-dimensional image in every area having a predeterminedsize.

As mentioned above, according to each Embodiment, it is possible tocalculate a displace amount by a least-squares method after eliminatinga solution which is intrinsically indeterminate. Moreover, according toone embodiment of the present invention, it is possible to correct analignment position between an optical image and a reference image byusing one of the more appropriate displacement amount in thedisplacement amount being the minimum SSD obtained as a result of theSSD calculation and the displacement amount obtained as a result of theleast-squares method calculation. That is, it is possible to correct analignment position by a more suitable optimization technique dependingupon an image. Thus, highly precise position alignment can be performed,thereby performing a high sensitive inspection.

What is expressed by “unit”, “circuit” or “step” in the abovedescription can be configured by a computer-executable program. They maybe executed by a software program or by any combination of software,hardware and/or firmware. Alternatively, they may be configured byhardware. When configured by a program, the program is recordable orstorable onto a recording medium, such as a magnetic disk drive,magnetic tape drive, FD or ROM (read-only memory). For example, thetable control circuit 114, the reference circuit 112, the comparingcircuit 108, the position alignment circuit 140, the respective circuitsin the position alignment circuit 140, and the like may be constitutedby electric circuits or the like. Alternatively, they may be achieved assoftware processable by the control computer 110, or achieved by acombination of electric circuits and software.

As mentioned above, Embodiments have been described with reference tothe concrete examples. However, the present invention is not limited tothese concrete examples. For example, transmitted light is used inEmbodiments, but reflected light may also be used, or transmitted lightand reflected light may be used simultaneously. The reference image isgenerated from design data, but alternatively, data of a same patternphoto-captured by a sensor such as a photodiode array may be employed.In other words, it is equally preferable to employ the die to dieinspection or the die to database inspection.

Moreover, though apparatus configurations, control methods, etc. whichare not directly required in explaining the present invention are notdescribed, a necessary apparatus configuration and a necessary controlmethod can be appropriately selected and used.

Furthermore, all of target workpiece inspection apparatuses, targetworkpiece inspection methods, image alignment methods, and positionaldisplacement estimating methods which have the elements of the presentinvention and which can be appropriately changed in design by a personskilled in the art are included in the spirit and scope of theinvention.

Additional advantages and modification will readily occur to thoseskilled in the art. Therefore, the invention in its broader aspects isnot limited to the specific details and representative embodiments shownand described herein. Accordingly, various modifications may be madewithout departing from the spirit or scope of the general inventiveconcept as defined by the appended claims and their equivalents.

1. A pattern inspection apparatus comprising: an optical image acquiringunit configured to acquire an optical image of pattern to be inspectedon which a pattern is formed; a reference image generating unitconfigured to generate a reference image to be compared with the opticalimage; a normal matrix element calculating unit configured to calculatea plurality of elements of a normal matrix for a least-squares methodfor calculating a displacement amount displaced from a preliminaryalignment position between the optical image and the reference image; apattern type estimating unit configured to estimate a type of a patternindicated by the reference image, by using some of the plurality ofelements of the normal matrix; a displacement calculating unitconfigured to calculate the displacement amount displaced from thepreliminary alignment position based on the least-squares method, byusing a normal matrix obtained by deleting predetermined elements fromthe plurality of elements of the normal matrix depending upon the typeof the pattern which has been estimated; a position correcting unitconfigured to correct an alignment position between the optical imageand the reference image to a position displaced from the preliminaryalignment position by the displacement amount; and a comparing unitconfigured to compare the optical image and the reference image whosealignment position has been corrected.
 2. The apparatus according toclaim 1, wherein the pattern type estimating unit estimates the type ofthe pattern indicated by the reference image by using a trace and anabsolute value of a determinant of a matrix of two rows and two columnscomposed of some of the plurality of elements of the normal matrix. 3.The apparatus according to claim 2, wherein as elements of the matrix oftwo rows and two columns, a total sum of first differentiation valuesobtained by space differentiating data of the reference image in Xdirection with respect to each pixel, a total sum of seconddifferentiation values obtained by space differentiating data of thereference image in Y direction with respect to each pixel, a total sumof squared first differentiation values, and a total sum of squaredsecond differentiation values are used.
 4. The apparatus according toclaim 3, wherein the pattern type estimating unit estimates the type ofthe pattern to be an no-structure pattern in which a whole image is apattern when the trace is smaller than a first threshold value and theabsolute value of the determinant is smaller than a second thresholdvalue, to be a line and space pattern with an arbitrary angle when thetrace is greater than or equal to the first threshold value and theabsolute value of the determinant is smaller than the second thresholdvalue, and to be other pattern when the trace is greater than or equalto the first threshold value and the absolute value of the determinantis greater than or equal to the second threshold value, and thedisplacement calculating unit, when the type of the pattern is estimatedto be the no-structure pattern, calculates the displacement amount byusing a normal matrix obtained by deleting elements using the firstdifferentiation value of the normal matrix and elements using the seconddifferentiation value of the normal matrix from the plurality ofelements of the normal matrix, and when the type of the pattern isestimated to be the line and space pattern, calculates the displacementamount by using a normal matrix obtained by deleting elements usingpredetermined one of the first differentiation value and the seconddifferentiation value of the normal matrix from the plurality ofelements of the normal matrix.
 5. The apparatus according to claim 1,wherein the displacement calculating unit divides the reference image byweighting and calculates the displacement amount for each of dividedareas.
 6. A pattern inspection apparatus comprising: an optical imageacquiring unit configured to acquire an optical image of pattern to beinspected on which a pattern is formed; a reference image generatingunit configured to generate a reference image to be compared with theoptical image; a first SSD (Sum of Squared Difference) calculating unitconfigured to calculate a first displacement amount from a preliminaryalignment position between the optical image and the reference image toa first position where an SSD between a pixel value of the optical imageand a pixel value of the reference image is minimized; a normal matrixelement calculating unit configured to calculate a plurality of elementsof a normal matrix for a least-squares method for calculating adisplacement amount displaced from the preliminary alignment position; apattern type estimating unit configured to estimate a type of a patternindicated by the reference image, by using some of the plurality ofelements of the normal matrix; a displacement calculating unitconfigured to calculate a second displacement amount displaced from thepreliminary alignment position based on the least-squares method, byusing a normal matrix obtained by deleting predetermined elements fromthe plurality of elements of the normal matrix depending upon the typeof the pattern which has been estimated; a second SSD calculating unitconfigured to calculate an SSD between the pixel value of the opticalimage and the pixel value of the reference image at a second positiondisplaced from the preliminary alignment position by the seconddisplacement amount; an SSD estimating unit configured to estimate whichof the SSD at the first position and the SSD at the second position issmaller; a position correcting unit configured to correct an alignmentposition between the optical image and the reference image to a positionwhere a smaller SSD as a result estimated by the SSD estimating unit isobtained; and a comparing unit configured to compare the optical imageand the reference image whose alignment position has been corrected. 7.The apparatus according to claim 6, wherein the pattern type estimatingunit estimates the type of the pattern indicated by the reference imageby using a trace and an absolute value of a determinant of a matrix oftwo rows and two columns composed of some of the plurality of elementsof the normal matrix.
 8. The apparatus according to claim 7, wherein aselements of the matrix of two rows and two columns, a total sum of firstdifferentiation values obtained by space differentiating data of thereference image in X direction with respect to each pixel, a total sumof second differentiation values obtained by space differentiating dataof the reference image in Y direction with respect to each pixel, atotal sum of squared first differentiation values, and a total sum ofsquared second differentiation values are used.
 9. The apparatusaccording to claim 8, wherein the pattern type estimating unit estimatesthe type of the pattern to be an no-structure pattern in which a wholeimage is a pattern when the trace is smaller than a first thresholdvalue and the absolute value of the determinant is smaller than a secondthreshold value, to be a line and space pattern with an arbitrary anglewhen the trace is greater than or equal to the first threshold value andthe absolute value of the determinant is smaller than the secondthreshold value, and to be other pattern when the trace is greater thanor equal to the first threshold value and the absolute value of thedeterminant is greater than or equal to the second threshold value, andthe displacement calculating unit, when the type of the pattern isestimated to be the no-structure pattern, calculates the seconddisplacement amount by using a normal matrix obtained by deletingelements using the first differentiation value of the normal matrix andelements using the second differentiation value of the normal matrixfrom the plurality of elements of the normal matrix, and when the typeof the pattern is estimated to be the line and space pattern, calculatesthe second displacement amount by using a normal matrix obtained bydeleting elements using predetermined one of the first differentiationvalue and the second differentiation value of the normal matrix from theplurality of elements of the normal matrix.
 10. The apparatus accordingto claim 6, wherein the first SSD calculating unit calculates the firstdisplacement amount by performing displacing by subpixel unit to aposition where the SSD is minimized.
 11. The apparatus according toclaim 6, wherein the displacement calculating unit divides the referenceimage by weighting and calculates the displacement amount for each ofdivided areas.
 12. A pattern inspection apparatus comprising: an opticalimage acquiring unit configured to acquire an optical image of patternto be inspected on which a pattern is formed; a reference imagegenerating unit configured to generate a reference image to be comparedwith the optical image; a first SSD (Sum of Squared Difference)calculating unit configured to calculate a first displacement amountfrom a preliminary alignment position between the optical image and thereference image to a first position where an SSD between a pixel valueof the optical image and a pixel value of the reference image isminimized; a normal matrix element calculating unit configured tocalculate a plurality of elements of a normal matrix for a least-squaresmethod for calculating a displacement amount displaced from the firstposition; a pattern type estimating unit configured to estimate a typeof a pattern indicated by the reference image, by using some of theplurality of elements of the normal matrix; a displacement calculatingunit configured to calculate a second displacement amount displaced fromthe first position based on the least-squares method, by using a normalmatrix obtained by deleting predetermined elements from the plurality ofelements of the normal matrix depending upon the type of the patternwhich has been estimated; a second SSD calculating unit configured tocalculate an SSD between the pixel value of the optical image and thepixel value of the reference image at a second position displaced fromthe first position by the second displacement amount; an SSD estimatingunit configured to estimate which of the SSD at the first position andthe SSD at the second position is smaller; a position correcting unitconfigured to correct an alignment position between the optical imageand the reference image to a position where a smaller SSD as a resultestimated by the SSD estimating unit is obtained; and a comparing unitconfigured to compare the optical image and the reference image whosealignment position has been corrected.
 13. The apparatus according toclaim 12, wherein the pattern type estimating unit estimates the type ofthe pattern indicated by the reference image by using a trace and anabsolute value of a determinant of a matrix of two rows and two columnscomposed of some of the plurality of elements of the normal matrix. 14.The apparatus according to claim 13, wherein as elements of the matrixof two rows and two columns, a total sum of first differentiation valuesobtained by space differentiating data of the reference image in Xdirection with respect to each pixel, a total sum of seconddifferentiation values obtained by space differentiating data of thereference image in Y direction with respect to each pixel, a total sumof squared first differentiation values, and a total sum of squaredsecond differentiation values are used.
 15. The apparatus according toclaim 14, wherein the pattern type estimating unit estimates the type ofthe pattern to be an no-structure pattern in which a whole image is apattern when the trace is smaller than a first threshold value and theabsolute value of the determinant is smaller than a second thresholdvalue, to be a line and space pattern with an arbitrary angle when thetrace is greater than or equal to the first threshold value and theabsolute value of the determinant is smaller than the second thresholdvalue, and to be other pattern when the trace is greater than or equalto the first threshold value and the absolute value of the determinantis greater than or equal to the second threshold value, and thedisplacement calculating unit, when the type of the pattern is estimatedto be the no-structure pattern, calculates the second displacementamount by using a normal matrix obtained by deleting elements using thefirst differentiation value of the normal matrix and elements using thesecond differentiation value of the normal matrix from the plurality ofelements of the normal matrix, and when the type of the pattern isestimated to be the line and space pattern, calculates the seconddisplacement amount by using a normal matrix obtained by deletingelements using predetermined one of the first differentiation value andthe second differentiation value of the normal matrix from the pluralityof elements of the normal matrix.
 16. The apparatus according to claim12, wherein the first SSD calculating unit calculates the firstdisplacement amount by performing displacing by subpixel unit to aposition where the SSD is minimized.
 17. The apparatus according toclaim 12, wherein the displacement calculating unit divides thereference image by weighting and calculates the displacement amount foreach of divided areas.
 18. An image alignment method for aligning anoptical image and a reference image for use in a comparing inspection ofpattern to be inspected on which a pattern is formed, the methodcomprising: calculating a first displacement amount from a preliminaryalignment position between the optical image and the reference image toa first position where an SSD (Sum of Squared Difference) between apixel value of the optical image and a pixel value of the referenceimage is minimized; calculating a plurality of elements of a normalmatrix for a least-squares method for calculating a displacement amountdisplaced from the first position; estimating a type of a patternindicated by the reference image, by using some of the plurality ofelements of the normal matrix; calculating a second displacement amountdisplaced from the first position based on the least-squares method, byusing a normal matrix obtained by deleting predetermined elements fromthe plurality of elements of the normal matrix depending upon the typeof the pattern which has been estimated; calculating an SSD between thepixel value of the optical image and the pixel value of the referenceimage at a second position displaced from the first position by thesecond displacement amount; estimating which of the SSD at the firstposition and the SSD at the second position is smaller; and correctingan alignment position between the optical image and the reference imageto a position where a smaller SSD as a result of the estimating of SSDis obtained, to output a result of the correcting.
 19. (canceled)
 20. Acomputer-readable recording medium with a program recorded for causing acomputer to execute processes of: storing an optical image and areference image for use in a comparing inspection of pattern to beinspected on which a pattern is formed, in a storage device; reading theoptical image and the reference image from the storage device, andcalculating a plurality of elements of a normal matrix for aleast-squares method, based on a preliminary alignment position betweenthe optical image and the reference image; estimating a type of apattern indicated by the reference image, by using some of the pluralityof elements of the normal matrix; and calculating a displacement amountdisplaced from the preliminary alignment position based on theleast-squares method, by using a normal matrix obtained by deletingpredetermined elements from the plurality of elements of the normalmatrix depending upon the type of the pattern which has been estimated,to output the displacement amount calculated.